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ABSTRACT 

We studied the effect of primordial non-Gaussianity with varied bispectrum shapes on the 
number counts of signal-to-noise peaks in wide field cosmic shear maps. The two cosmolog- 
ical contributions to this particular weak lensing statistic, namely the chance projection of 
Large Scale Structure and the occurrence of real, cluster-sized dark matter halos, have been 
modeled semi-analytically, thus allowing to easily introduce the effect of non-Gaussian initial 
conditions. We performed a Fisher matrix analysis by taking into account the full covari- 
ance of the peak counts in order to forecast the joint constraints on the level of primordial 
non-Gaussianity and the amplitude of the matter power spectrum that are expected by future 
wide field imaging surveys. We find that positive-skewed non-Gaussianity increases the num- 
ber counts of cosmic shear peaks, more so at high signal-to-noise values, where the signal is 
mostly dominated by massive clusters as expected. The increment is at the level of ~ 1 % for 
/nl = 10 and ~ 10% for /nl = 100 for a local shape of the primordial bispectrum, while 
different bispectrum shapes give generically a smaller effect. For a future survey on the model 
of the proposed ESA space mission Euclid and by avoiding the strong assumption of being 
capable to distinguish the weak lensing signal of galaxy clusters from chance projection of 
Large Scale Structures we forecasted a 1 -cr error on the level of non-Gaussianity of ~ 30-40 
for the local and equilateral models, and of ~ 100 - 200 for the less explored enfolded and 
orthogonal bispectrum shapes. 

Key words: cosmology: theory - gravitational lensing: weak - cosmological parameters - 
large-scale structure of the Universe 



1 INTRODUCTION 

Gravitational lensing is one of the most powerful means of astro- 
physical investigation. While baryonic tracers of the dark matter 
distribution commonly rely on strong simplifying assumptions, the 
gravitational deflection of light tBarte l man n & Schneid er 2001) is 
sensitive solely to the overall incidence of any matter component 
along the line of sight. In practical applications, the occurrence of 
strongly distorted images in the core of massive galaxy clusters can 
return important information about the inner structure of dark mat- 
ter halos, while the weak systematic image distortion of large num- 
bers of background galaxies allows to efficiently trace the outskirts 
of clusters and the Large Scale Structure (LSS henceforth) in gen- 
eral. 

Cosmic shear, that is weak gravitational lensing on cosmo- 
logical scales, measures the redshift evolution of the LSS in the 
Universe weighted by a specific combination of angular diameter 
distances. Hence it combines cosmological tests based on the ge- 
ometry of the Universe with tests based on the growth of structures. 



Therefore, cosmic shear has long been recognized as an important 
tool for cosmology. In this paper we considered one particular cos- 
mic shear statistic, namely the abundance of signal-to-noise (S/N 
henceforth) ratio peaks in wide field weak lensing maps. Besides 
the intrinsic ellipticity distribution of background sources, peak 
counts are determined by the cosmological information encapsu- 
lated by the cosmic matter density fluctuations from linear to non- 
linear scales, namely LSS fluctuations and cluster-sized dark matter 
halos. Counting cosmic shear peaks without addressing their origin 
is not only easier to do than, for instance, counting weak lensing 
selected galaxy clusters, but also contains more cosmological lever- 
age (Dietrich & Hartlap 2009). 

In the present work we applied the abundance of S/N peaks in 
cosmic shear maps to the issue of cosmological initial conditions, 
and in particular whether or not the primordial (dark) matter density 
fluctuations are distributed according to a Gaussian distribution. It 
is now well established that primordial non-Gaussianity has a con- 
siderable impact on different aspects of structure formation. Specif- 
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ically, a positively (negatively) skewed distribution of initial mat- 
ter fluctuations would produce both a more (less) abundant galaxy 
cluster population and more (less) biased structures with respect to 
the underlying matter density field. Therefore, it is expected that 
the peak statistics would return valuable constraints on the level 

and shape of primordial non-Gaussianity. 

Besides several pioneering works jMessina et iD Il99d : 
iMoscardini et al.lll99ll : [Weinberg & C oi3 ll992|) , the problem of 
constraining deviations from primordial Gaussianity by means 
different from the Cosmic Microwave Background (CMB) 
intrinsic anisotropies has recently attracted renewed attention 
in the literature, with efforts dir ected towards the abundance 
of non-lin e ar str u ctures cMataiTese. Verde. & Jimene j [2OO0I: 



Verde et alJ I2OO0I; iMathis Diego, & Silkl |2004 iGrossi et al 



2OO71, 12009; 'Maggiore &. Riotto"2010c'), halo biasing /Dalai et al.l 
2008; McDonald 2008; Fedeli. Moscardini. & Matarrese 2009; 
Fedeli et al.l 12 010). galaxy bispectrum dSefusatti & Komatsul 
20071; IJeong & Komatsu .2009,) , mass density distribution 



jGrossi et aljboOSi) and topology iMatsubaral 20031 ; iHikage et al.l 
l2008b . cosmic shear l lFedeli & Moscar dini' '201 Ol, Pace et al . 
2010), integrated Sachs-Wolfe eff ecM|Afshordi & TollevI |2008| ; 
ICarbone. Verde. & MatarreselbOOSi), Lya flux from low-de nsity in- 
tergalactic medium I'Viel et al. "200% 21 — cm fl uctuations jCooravl 
[2OO6; Pillepich, Porciani, & Matar res3 I2OO7I) and reionization 
l lCrociani et al.ll2009i) . 

As a specific application of our investigation on cosmic shear, 
we refer to a future half-sky optical/near-infrared imagin g survey 
on th e model of the ESA Cosmic Vision proposal Euclid 1 LaureiisI 
I2OO9I) . Our aim is to forecast the constraints on the level of pri- 
mordial non-Gaussianity that are expected by counting S/N peaks 
in cosmic shear maps that should be produced by Euclid. We shall 
also investigate how these constraints change upon modification of 
several survey parameters and detection criteria, such as the imag- 
ing depth, the scale of optimal filtering, and the S/N threshold. 

Throughout this work we adopted for the reference Gaussian 
model the cosmolog ical parameter set gi ven by the latest analysis 
of the WMAP data jKomatsu et ai]|201lh . namely the matter den- 
sity parameter n„.o = 0.272, the cosmological constant density pa- 
rameter JIa.o = 0.728, the baryon density parameter fib.o = 0.046, 
the Hubble constant h = //o/(100kms"' Mpc"') = 0.704, and the 
matter power spectrum normalization set by erg = 0.809. The rest 
of our paper is organized as follows. In Section |2] we summarize 
the various non-Gaussian cosmologies adopted, while in Section[3] 
we recall the modifications that these introduce to the cluster mass 
function, the large scale bias, and the matter power spectrum. In 
Section|4]we describe the various contributions to the weak lensing 
peak number counts and how they have been modeled. In Section|5] 
we show the effect of primordial non-Gaussianity on the abundance 
of cosmic shear peaks, while in Section|6]we report the Fisher ma- 
trix analysis that we performed in order to forecast constraints on 
the initial conditions. Finally in Section|7]we summarize our con- 
clusions. 



2 NON-GAUSSIAN COSMOLOGIES 



Exten s ions of the most standard model of inflation jStarobinskiil 
1 19791 ; iGutlj Il98ll ; iLindd Il982h can produce substantial de- 
viations from a Gaussian di stribution of primo r dial density 
and potential fluctuations (see iBartolo et |2004| ; IChenI I2OI0I ; 
iDesiacgues & Seliakll2010l for recent reviews). The amount and 
shape of these deviations depend critically on the kind of non- 



standard inflationary model that one has in mind, as will be detailed 
later on. 

A particularly convenient (although not unique) way to de- 
scribe generic deviations from a Gaussian distribution consists in 
writing the gauge-invariant Bardeen's potential O as the sum of a 
Gaus s ian random field and a quadratic correct ion ( Salopek & Bond 
1990| ; iGanguiet alii 1994 I Verde et al.1 [2OO0I ; Ixomatsu & Spergel 
200 ih . according to 

t = <I>G + /nL * (Og - <*G>) ■ (1) 

On sub-horizon scales the Bardeen's potential equals minus the 
Newtonian peculiar potential. The parameter /nl in Eq. (TJ deter- 
mines the amplitude of non-Gaussianity, and it is in general de- 
pendent on the scale. The symbol * denotes convolution between 
functions, and reduces to standard multiplication upon constancy 
of /nl. In the following we adopted the large-scale structur e con- 
vention (as opposed to the CMB convention, see[A fshordi & TollevI 
20081; ICarbone et al.ll2008l ; IPillepich. Porciani. & Hahn 20091 and 
Grossi et al .I2OO9I) for defining the fundamental parameter /nl. Ac- 
cording to this, the primordial value of O has to be linearly extrap- 
olated at z = 0, and as a consequence the constraints given on /ml 
by the CMB have to be rais ed by ~ 30 per cent to comply with this 
paper's convention (see also lFedeli. Moscardini. & Matarresel2009l 
for a concise explanation). 

In the case in which /nl t the potential O is a random field 
with a non-Gaussian probability distribution. Therefore, the field 
itself cannot be described by the power spectrum Pdk) = Bld^^* 
alone, rather higher-order moments are needed, for instance the bis- 
pectrum Bf,(k[ , k2, A3). The bispectrum is the Fourier transform of 
the three-point correlation function {(t>{ki)(I>(ki)<t>(ki)) and it can 
hence be implicitly defined as 

{(S>(k,mki)(S>{ki)) = (2nf6D{k,+k2 + ki)B^{k,,k2,ki) . (2) 

As mentioned above understanding the shape of non- 
Gaussianity is of fundamental importance in order to pinpoint the 
physics of the early Universe and the evolution of the inflaton field 
in particular. For this reason, in this work we considered four differ- 
ent shapes of the potential bispectrum, arising from different modi- 
fications of the standard inflationary scenario. We summarize them 
in the following, refe rring the reader to the quoted references and 
to lFedeli et ^ hoiOh for further details. 



Local shape 

The standard single-field inflationary scenario generates negligibly 
small deviations from Gaussianity. These deviations are said to be 
of the local shape, and the related bispectrum of the Bardeen's po- 
tential is maximized for squeezed configurations, where one of the 
three wavevectors has much smaller magnitude than the other two. 
In this case the parameter /nl is usually assumed to be a constant, 
and it is expected to be of the same order of the slow-roll param- 
eters (Falk, Rangarajan, & Srednicki 1993), that are very close to 
zero. 

However non-Gaussianities of the local shape can also be gen- 
erated in the case in which an additional light scalar field, differ- 
ent from the inflaton, contributes to the observed curvature pertur- 
bations tBabich. Crem inelli. & Zaldarriaga 2004). This happens , 
for instance, in curvaton models jSasaki. Valiviita. & Wand sliooel ; 



[Assadullahi, Valiviita, & Wandsll2007h or in multi-fields models 
(iBartolo, Matarrese, & Riottd l2002l ; lBemardeau & Uzan|[20o3) . In 
this case the parameter /nl is allowed to be substantially different 
from zero. 
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Equilateral shape 

In some inflationary models the kinetic term of the inflaton 
Lagrangian is not standard, containing higher-order derivatives 
of the field itself. One signifi cant example of this is the DBI 
mode l ( lAlishahiha. Silverstein. & To ng 2004; iSilverstein & Tond 
2004 see also lArkani-Hamed et all i2004: .Seerv & LidsevUlOOsT 



Li. Wang. & Wandl2008h . In this case the primordial bispectrum 



is maximized for configurations where the three wavevectors have 
approximately the sam e amplitude, and i t is we ll represented by the 
template introduced bv lCreminelli et alJ (l2007h . 

Given this, we have the freedom to insert a running 
y(ki,k2,kj) for /nli since this parameter is not forced to be 
constant in this case. The form we chose reads dCherj I2OO5I : 
ILo Verde et al.ll2008l : ICrociani et al.ll2009i) 



(3) 



In all calculations that follow we considered this running as part of 
the equilateral bispectrum. This is important since different authors 
choose different runnings, or no running at all, for the equilateral 
shape. We adopted the exponent k = -0.2 in the remainder of this 
work, that increase the level of non-Gaussianity at scales smaller 
than that corresponding to kcuB = 0.086/? Mpc"'. This coincides 
with t he larger multi pole used in the CMB analysis by the WMAP 
team jKomatsu et ^,.2009 . 201 1.^ . ( ~ 700. 



Enfolded shape 

For deviations from Gaussianity evaluated in the regular Bunch- 
Davies vacuum state, the primordial potential bispectrum is 
of local or equilateral shape, depending on whether or not 
higher-order derivatives play a significant role in the evolu- 
tion of the inflaton field. If the Bunch-Davies vacuum hy- 
pothesis is dropped, the resulting bi s pectru m is maximal for 
squas hed configurations (Chenetal. 2007; Holman & TollevI 
|2008|) . iMeerburg, van der Schaar, & Corasaniti (2009) found a 
template that describes very well the properties of t his enfolded- 
shape bispectrum (see however fCreminelli et al.120101 for a slightly 
different template of the physical model). 



Orthogonal shape 

A shape of the bispectrum can be constructed that is 
nearly orthogonal to both the lo cal and equilateral forms 
dSenatore. Smith. & Zaldarriagall2OI0l) . Constraints on the level of 
non-Gaussianity compatible with the CMB in the local, equilateral 
and o rthogonal scenarios were recently given by the WMAP 
team teomatsuet"ai]|20Ilh . while constraints on enfolded non- 
Gaussianity from galaxy bias were given by IVerde & Matarres^ 
( l2009l) 

Although there is no theoretical prescription against a 
rurming of the /nl parameter with the scale in the enfolded and 
orthogonal shapes, we decided not to include one. The reason 
for this is that there is no first principle that can guide one in 
the choice of a particular kind of running, and until now no 
work has addressed th e proble m of a running for these shapes 
jpergusson & Shellard |2009| ; iFergusson. Liguori. & Shellar j 
l2010 f). Moreover, we recall that the four non-Gaussian shapes 
described above are not independent, rather the orthogonal shape 



can be obtained as a suita ble linear combination of the equilat- 
eral and enfolde d shapes ( Senatore, Smith. &Zaldarriag3l20ld; 
Wagn er. Verde. & Boubekeur 2010). Nevertheless, we performed 
computations for the orthogonal model as well, since it gives a 
different signature on the evolution of the LSS. 



3 COSMOLOGICAL OBSERVABLES 

Primordial non-Gaussianity produces modifications in the statistics 
of density peaks, resulting in differences in the mass function of 
cosmic objects and the bias of dark matter halos with respect to the 
underlying smooth density field. In the following we summarize 
how these modifications have been taken into account in the present 
work, and how they reflect into modifications to the matter power 
spectrum. 



3.1 Mass function 

For the non-Gaussian modification to the mass function of cos - 
mic objects we adopted the prescription of ILo Verde et alj ( l2008h . 
The main assumption behind it is that the effect of primor- 
dial non-Gaussianity on the mass function is independent of the 
prescription adopted to describe the mass function itself. This 
means that, if rt<,°'(M,z) and «ps(M,z) are the non-Gaussian and 
Gaussian mass functions , respectively, computed according to the 
iPress & Schechted ( [T974h formula, we can define a correction fac- 



tor ??(M,z) = nf^{M,z)lnps{M,z). Then, the non-Gaussian mass 
function computed according to an arbitrary prescription, n{M, z) 
can be related to its Gaussian counterpart through 

n{M,z) = 'R{M,z)n'^°\M,z} 



In orde r 
ILo Verde et al 



(4) 

to compute «ps(M, z), and hence 'R(M,z), 

1 20081) perfo rmed an Edgeworth expansion 



dSlinnikov & MoessneJ Il998l) of the probability distribution for 
the smoothed density fluctu ations field, tru ncating it at the linear 
term in ctm- The resulting IPress & Schechtej jl974h -like mass 
function reads 



;jps(M,z) 



/f Pm 
S n M 



exp 



6 

1 dSj 
6 dM 



St(z) 
^4 



Sl(z) 



d\wcr„ (6^{z) 



dM 



1 + 



Sliz) 



(5) 



In the previous equation p,^ = 3HgQ.^Q/SnG is the comoving 
matter density in the Universe, ctm is the rms of density fluc- 
tuations smoothed on a scale corresponding to the mass M, and 
S^iz) = A(./D+(z). The function S^IM) = iijiM)/^*^ is the reduced 
skewness of the non-Gaussian distribution, and the skewness ^ii(M) 
can be computed as 



y"3(M) 



f. 



M«(^,)M«(fc)M«fc)x 



X <Cl)(ki)0(k2)0(k3))- 



<ikidk2<ik: 



(6) 



The last thing that remains to be defined is the function MR{k), that 
relates the density fluctuations smoothed on some scale R to the 
respective peculiar potential, 

^ T{k)k^ 



MK{k) : 



3 //gfin 



-WRik); 



(7) 
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where T{K) is the matter transfer function and W[((k) is the Fourier 

transform of the top-hat window function^ 

In this work we adopted the Bar deen et al matter 
transfe r function, with the shape factor correction of Sugivamj 
( Il995t). This reproduces fa irly well the more sophisticated recipe 
of lEisenstein & Hul ( Il998h except for the presence of the baryon 
acoustic oscillation, that anyway is not of interest here. We ad- 
dit ionally adopted as refe rence mass function th e prescription 
of ISheth & TomienI f2Cio3) (see I Jenkins et al.ll200ll : [Warren et all 
I2OO6I : Tinker et"H]|2008l for alternative prescriptions). Other ap- 
proaches also exist for computing the non-Gaussian correction to 
the mass funct ion, that give restilts in broad agreem ent with those 
obtained here l lMatarrese. Verde. & Jimeneal2000h . In computing 
the non-Gaussian corrections to the mass function we have taken 
into account t he correction t o the critical overdensity for collapse 
suggested by iGrossi et"aL I ( l2009h (see also Maggiore & Riottg 
according to which — > -y/^, with q ~ 0.8. 



3.2 Halo bias 

Recently much attention has been devoted to the effect of 
primordial non-Gaus sianity on halo bias, and the use thereof 
for c onstraining /nl (Matarrese & Verde 2008; Dalai et al."200d; 
IVerde & Matarrese 2009; Carbone, Verde, & Matarrese 2008). In 
particular, these works have shown that primordial non-Gaussianity 
introduces a scale dependence on the large scale halo bias. This pe- 
culia rity allows to place already stringent constr aints from existing 
data jSlosar et al.ll2008l ; lAfshordi & ToIIevll2008 ). 

The non-Gaussian halo bias can be written in a relatively 
straightforward way in terms of its Gaussian counterpart as 
dCarbone. Mena. & Verd3l2010l) 

KM, z, k) = fo<°'(M, z) +l3R(k)crl [b^°\M, z) - if , (8) 

where the function f^nik) encapsulates all the scale dependence of 
the non-Gaussian correction to the bias, and reads 



8;rV^M 



(9) 



where a = k^ + + 2k(fi. In the simple case of local bispectrum 
shape it can be shown that the function Pit(k) should scale as oc k^- 
at large scales, so that a substantial boost (if /nl > 0) in the halo 
bias is expected at those scal es. For the Gaussian bi a s W^\ M, z) we 
adopted the prescription of ISheth. Mo. & TormenI j200lh . In this 
case, since the correction to the Gaussian bias is written in term 
of the Ga ussian bias itself, th e ellipsoidal collapse correction sug- 
gested bv lGrossi et ^ {2009) is not necessary. 

It is interesting to note that, while for the first three non- 
Gaussian shapes introduced in Section (2] a positive /nl implies 
both a positive skewness of the matter density field and a positive 
correction to the large scale halo bias, the opposite is true for the 
fourth shape, the orthogonal one. This is a fact to be kept in mind 
when interpreting our results in the subsequent Sections. 



3.3 Matter power spectrum 

Another cosmological observable that is relevant for our purposes 
and gets modified in case of primordial non-Gaussianity is the mat- 
ter power spectrum. The latter requires a little bit more of care 



with respect to the mass function and linear bias. Particularly, we 
should have a reliable way to parametrize the fully non-linear three- 
dimensional power spectrum, since part of the lensing signal is 
given by the integral thereof along t he line of sight. In order to 
do that we followed the approach of iFedeli & Moscardinil j2010l) 
(and references therein) and made use of the halo model in order to 

represent the matter power s pectrum . 

The halo model jSeliakI I2OO0I ; iMa & Frvl I2OO0I ; 
ICoorav & ShethI l2002h is a physically motivated framework 
that allows to compute the correlation function of various LSS 
tracers, including the dark matter particles themselves. It is based 
on the simple consideration that if all the matter in the Universe is 
locked into halos, then the contribution to the correlation function 
comes from particle pairs sitting in separated halos at large scales 
and from particle pairs belonging to the same halo at small scales. 
Accordingly, the power spectrum can be written as the sum of two 
terms describing the two contributions, P{k, z) = Pi(k, z) + Pi{k, z), 
with 



Pi{k,z) 



and 



PAKz) 



: r ;j(M, 
Jo 



z) 



p(k,M,z) 



dM . 



r°° Pik, M, z) 

n(M,z)b(M,z,k)— — —dM 
Jo Pm 



PL(k,z). 



(10) 



(11) 



In the previous set of equations P]^{k, z) is the linear matter power 
spectrum and p(k, M, z) is the Fourier transform of the mean dark 
matter halo density profile which is also in principle modified 
by pri mor dial non-Gaussianity, as sugge s ted in lAvila-Reese et al] 
( I2OO3I) and lSmith. Desiacques. & MarianI j2010|). Ne v erthel ess, we 
adopted a standard iNavarro. Frenk. & White! (ll996L Il997h shape 
(NFW henceforth) for p(k, M, z) because the functional forms for 
the modified density profiles are proven to hold only for the local 
shape, and it is not clear whether this is the case for other shapes 
as well. Therefore we ignored these changes, making our analysis 
more conservative, since a positive (negative) /nl implies a more 
(less) centrally concentrated dark matter profile, and consequently 
an increment (decrement) in the number of peaks. 

For practical details on the i mplementation of the hal o 
mode l we refer the intereste d reader to lAmara" & Refregierl ( l2004h ; 
Fedel i & Moscardinil ilOUt) . Here we just note that the mass func- 
tion and the halo bias, both entering in the halo model, are modified 
by primordial non-Gaussianity according to Sections [3 . 1 1 and [T2l 



4 MATTER FLUCTUATIONS WITH COSMIC SHEAR 

In this paper we focused on the number counts of weak lensing 
peaks, that directly reflect the distribution of dark matter fluctua- 
tions. Counting directly the peaks allows to minimize the physi- 
cal assumptions necessary in going from the adopted cosmological 
model to the data outcome. As a matter of fact we based the final 
results on dark matter physics only, by observing quantities insensi- 
tive to baryon physics and by relating the model prediction directly 
to the measured S/N ratios. In particular, the latter point avoids the 
difficult task of disentangling in a clear way the intrinsic nature of 
each single peak: noise fluctuation, LSS line-of-sight superimposi- 
tion or actual galaxy cluster. 
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4.1 Lensing formalism 

The deflection properties of isolated lenses are fully described by 
their two-dimensional lensing potential, 

2 



2 flds r 

c- DiD, Jo 



<b{Die,z)dz ., 



(12) 



where O is the Newtonian gravitational potential and D^, and 
Dds are the angular-diameter distances between the observer and 
the source, the observer and the lens, and the lens and the source, 
respectively. Finally, i- represents the physical distance out to the 
source sphere. 

The potential i/^ relates the angular positions j8 of a source and 
6 of its images on the observer's sky through the lens equation, 
fi = - Vtl/. For sources such as distant background galaxies it is 
possible to linearize the lens equation such that the induced image 
distortion is expressed by the Jacobian matrix 



(1 - k) 



1-gi 
-82 



-gl 
1 +gl 



(13) 



In the previous equation, k = is the convergence, responsi- 

ble for the isotropic magnification of an image relative to its source, 
and g = y/(l - k) is the reduced shear, quantifying the observed 
distortion. Here, 71 = {if/ n - if/ 22) /2 and 72 = 't',12 are the two 
components of the complex shear, where a comma denotes stan- 
dard differentiation with respect to coordinates on the lens plane. 
It is important to recall that since the angular size of the sources 
is unknown, only the reduced shear can be estimated starting from 
the observed ellipticity of the images. 

4.2 Contributions to the lensing signal 

The abundance of S/N peaks in cosmic shear maps is determined 
by the sum of three separate contributions, that we describe in the 
following. 

(i) The signal due to the occurrence of real non linear struc- 
tures such as galaxy clusters. Their abundance is evaluated in Sec- 
tion | 3.1| and their shear profil e can be ev aluated analytically (see 
iBarte lmann 1996; Meneghetti et alJI 20031) . Since the typical sepa- 
ration between the observed background galaxies is larger than the 
typical radius of the galaxy clusters' critical curves we can assume 
y - g throughout. 

(ii) The lensing signal due to chance projections of the LSS, 
given by the effective convergence power spectrum 

PAi) = — — J dw ^ Pi ^r^,w| , (14) 



where P(k,z) is the fully non-linear matter power spectrum dis- 
cussed in Section [373] £l,n_o is the present-day matter-density pa- 
rameter, Hq is the Hubble constant, c is the speed of light, a(w) is 
the scale factor at the comoving distance w, wh is the comoving dis- 
tance to the horizon, is the comoving angular-diameter distance, 
and W(w) is the weight function encorporating the line-of-sight in- 
tegral over the distribution of background sources G(w), 



W{w) - 



I dw G{w )- 



(15) 



/if(w') 

In our application we used a single (the tangential) component of 
the shear, for which the power spectrum reads Py^(£) = Pail)!!. 

(iii) The observational noise contribution from the intrinsic el- 
lipticity and finite number of background galaxies used to measure 
the shear signal, which has the white noise power spectrum 



1 

P.{e) = ~ , 

I n„ 



(16) 



determined by the number density Wg of galaxies suitable for weak 
lensing measurements and the variance o"| of their intrinsic ellip- 
ticity distribution. 

Since contributions (ii) and (iii) are well represented by two 
independent Gaussian random fields, their total power spectrum, is 
the sum of the two contributions, P(() = Py^(() + P^f). 



4.3 Optimal lensing filtering 

A method to measure gravitational lensing signatures in 
optical/near-infrared c atalogues of galaxies is the aperture mass 
l lSchneideret"ai] |l998l). It is a weighted average of the tangential 
component of the shear 7, over the position on the sky. 



A(0) 



^ r d'0'rt(0'. 



e)Q(\\0' - 0\\) , 



(17) 



where the radial filter function Q determines the statistical 
properties of the estimate A we are interested in. In this paper 
we aim at constraining the deviations from primordial Gaussianity, 
which mostly affect the non-linear part of structure formation and 
thus the high mass end of the mass function. Therefore , we adopted 
the linear matched filter defined bv lMaturi et al J ( [20051) which is de- 
signed to maximize the S/N ratio of non linear structures such as 
galaxy clusters, i.e those structures which are most sensitive to /nl. 



Qit) = a 



m 



with a = 



r 



d'e 



(18) 



Pit) Jk2 P{{) 

Here the filter is represented in the Fourier space for convenience, 
f{{) is the Fourier transform of the expected shear profile of indi- 
vidual halos, in our case the Fourier transform of the NFW shear 
profile, and P{i) = Py,{{) + Ps(.0 described in Section|42] Note that 
the shape of the optimal filter depends on the assumed parameters 
for the NFW density profile, particularly on the angular size corre- 
sponding to the scale radius 9^. Moreover, since the filter is only a 
radial function, its Fourier transform depends only on£ = \\t\\ 

We further define the variance of the aperture mass estimate A 

as 



r+°° £d{ , , 



(19) 



so that it is related only to the non-cosmological signal component, 
i.e. Pf. It is worth mentioning that in practical applications A is 
approximated as 



1 " 

A{0) = - V eu(0)Qm ■ 
n 



0\\) ■ 



(20) 



where ei,,(0) is the tangential ellipticity of a galaxy image located 
at the position S, with respect to 0, and which provides an estimate 
for 7t. 



4.4 Expected shear peaks number counts 

Once the lensing signals and the observational strategy have been 
defined (see Sections [4.2l and l431 l. the weak lensing number counts 
can be predicted as the sum of two independent components. On 
one hand the non-linear structures, i.e. galaxy clusters, whose de- 
tection depends on their intrinsic abundance given by the mass 
function discussed in Section [STT] and by their observed S/N ra- 
tio depending on their expected signal A (see Eq.ll7t and variance 
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Figure 1. Left panel. The different contributions to the abundance of shear S/N peaks in the reference Gaussian model. The black short dashed line shows the 
contribution given by LSS projection and data noise, while the blue long dashed line refers to the contiibution of the real dark matter clumps. The heavy solid 
red line is the total count distribution. We adopted optimal filtering with a scale radius of 6s = 2' . Right panel. The total number of shear S/N peaks for three 
different filter scale radii, as labeled. 



cr\ (see Eq. |19t . On the other hand the lensing LSS and instm- 
mental noise which are well approximated by a Gaussian random 
field and for which their number counts above a specific threshold 
S/N,i, = jti, can be easily predicted, 

= f ^ exp (- A) , (21) 

as shown bv lMaturi et alj ( [201(]h . Here cti, defined as 

Jo — ' (22) 

is the LSS plus noise lensing field variance which depends on 
the observational noise, the convergence power spectrum and the 
adopted filter. Note that, by approximating this therm as a Gaus- 
sian random field, the primordial non-Gaussianity is accounted 
through the LSS power spectrum only so that its leverage is not 
fully included. For our purpose this is a negligible approxima- 
tion since we anyway expect this effect to be small and put our 
final predictions on a conservative side. These two contributions 
carry complementary cosmological information, thus the statistics 
of cosmic shear peaks are in principle valuable cosmological tools 
jPietrich & HartladI 20091) . 

The two cosmological contributions to the abundance of cos- 
mic shear peaks are exemplified in the left panel of Figure[T]for the 
reference Gaussian cosmological model. For this Figure and in the 
rest of the work we ado pted the sourc e redsh ift distribution obser- 
vationally derived by iBeniamin et al.l ( l2007h . while their average 
surface number density is Wg = 40 arcmin"^, a level achieved by 
actual weak lensing surveys and which should be reached by EU- 
CLID (a prediction fo r a sample of exist ing and future weak lens- 
ing surveys is given bv lMaturi et alj2010l) . As expected, the chance 
alignment of the LSS dominates the counts at small S/N ratio val- 
ues, while the contribution coming from the occurrence of real dark 



matter clumps starts to be relevant at around S/N~ 4 and dominates 
at higher S/N values. We additionally show in the right panel of the 
same Figure the total number counts expected for three different 
scale radii of the optimal filter. As it can be noted, the larger the 
scale, the larger the contribution of galaxy clusters with respect to 
the LSS and the noise, although the latter remains always dominant 
at smaller, but still large, S/N. This behavior is expected because 
for larger filter scales the detections will have larger areas, thus in- 
creasing their blending where their number density is higher, typi- 
cally at lower S/N. Also, the overall S/N ratio values at higher S/N 
are expected to be larger, since a larger filter scale implies a better 
statistic for the background galaxies reducing their shot noise and 
intrinsic ellipticity. These features are ensured by the adopted op- 
timal filter which maximizes the cluster signal against the one of 
LSS and noise. In the remainder of this paper we shall refer mainly 
to a scale radius of 2', unless explicitly noted otherwise. 



5 SIGNATURE OF PRIMORDIAL NON-GAUSSIANITY 

The effect of primordial non-Gaussianity on the statistics of shear 
peaks is exemplified in Figure|2]for the four different non-Gaussian 
shapes that were discussed in Section |2l and for levels of non- 
Gaussianity /nl = ±50 and /nl = ±100. As it can be seen, the 
effect of non-Gaussianity is more evident at high values of S/N ra- 
tio, where it is mostly given by the effect on the mass function of 
dark matter halos. We remind the reader that our predictions re- 
garding the impact of primordial non-Gaussianity at the S/N ratios 
dominated by the LSS are conservative as discussed in Section |4!4l 
Also, as expected, the effect is largest for local non-Gaussianity for 
which, at high S/N, the counts can be modified by up to ~ 7% for 
/nl = ±100 and smallest for the orthogonal shape, where it barely 
reaches ~ 1%. For the orthogonal shape positive values of /nl pro- 
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Figure 2. The ratio of the peak number counts in different models with primordial non-Gaussianity to the same quantities in the reference Gaussian ACDM 
cosmology, as a function of the S/N ratio. Four different non-Gaussian shapes and four different values of /nl are shown, as labeled in the plot. 



vide a decrement in the number counts of shear peaks at high S/N 
values, while the opposite is true for all other models. This agrees 
with the behavior of the skewness that has been discussed in Sec- 
tion [2] Although the effect of non-Gaussianity is substantially re- 
duced for shapes different from the local one, it is likely that use- 
ful constraints can be put on /nl for these shapes as well. For in- 
stance in the orthogonal case, is constrained only a t the level 
of a few hundreds by CMB data jKomatsu et al. 2011 ). O ur pre - 
diction is compatible with the recent work of iMarian et all jZOlOl) . 
who found an effect on the shear peaks number counts at the level 
of ~ 10% for the highest significance detections. Although a direct 
comparison of the two works is not possible due to different source 
redshift distributions, different detection schemes, etc., we actually 
verified that by reducing the scale radius of the optimal filter down 
to 6s = r the effect is raised to ~ 10%. Thus, we conclude that the 
two works are in broad agreement for the local non-Gaussian model 
concerning the magnitude of the effect. This is reassuring, since 
the two results are based on quite different premises and adopt dif- 
ferent methodologie s: we opted for semi-analytic modeling, while 
iMarian et all ( |2010|) resorted to the use of cosmological simula- 
tions. 



In Figure [3] we show the ratio of the peak number counts in 
non-Gaussian cosmologies with local shape (/nl = ±10) to the 
same quantities in the reference Gaussian case adopting three dif- 
ferent values of the power spectrum amplitude erg for all models. 
We considered ~ 2% fluctuations in the value of erg, comparable 
with the 1 - tr uncertainty around the fiducial value derived by the 
WMAP-7 release. Quite unexpectedly, the effect of primordial non- 
Gaussianity is higher for lower values of erg- We interpret this fact 
as due to the delayed structure formation implied by a lower nor- 
malization. This allows the effect of primordial non-Gaussianity 
on the matter density field to be retained for a longer time by the 
LSS, with the consequence that the S/N peak number counts are 
slightly more affected. In other words, since with a lower erg the on- 
set of non-linear evolution is delayed, gravitational clustering has 
less time to erase the cosmological initial conditions. 
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Figure 3. The ratio of the peak number counts in two different non- 
Gaussian models with local bispectmm shape to the same quantities in 
the reference Gaussian ACDM cosmology, as a function of the S/N ratio. 
The upper set of curves refer to /nl = 10, while the lower set refers to 
/nl = ~10. Different colors refer to different values of the matter power 
spectrum normalization erg, as labeled. 



6 COSMOLOGICAL CONSTRAINTS 

In this section we use a Fisher-matrix analysis to estimate the confi- 
dence level in the erg - /nl plane, as obtained with the weak lensing 
number counts. 



6.1 Fisher matrix analysis 

The Fisher matrix is given by 



(23) 



where X is the logarithm of the likelihood function and ^ = 
(o"8i,/nl) are the free parameters of the number counts model. In 
case of a multivariate Gaussian likelihood, the Fisher matrix can be 
written as 



F,j = ^Tr[A,Aj + C-'M,j] , 



(24) 



where A; = C^'Cj, M,, = 2(dfi/d^d(did/d^j), C is the data co- 
variance matrix, and yu is the assumed model, i.e. the peaks number 
counts for different S/N. Here a comma denotes differentiation with 
respect to the relevant cosmological parameter. Since C has a very 
weak dependence on the model parameters erg and /nl the first term 
in Eq.|24]is negligible because of C,,. We evaluated the Fisher ma- 
trix at the fiducial point erg = 0.809 an d /nl = as measur ed by 
WMAP-7 for a standard ACDM model feomatsu et^l201lh . 

As discussed in Section |4!4l the number counts model fi is 
given by two components. On one hand, the non linear structures 
whose number counts are independent for different S/N ratio bins 
and whose covariance contribution is therefore diagonal with a 
Poissonian amplitude. On the other hand, the contribution given 
by the LSS and noise for which the number counts are relevant for 




Figure 4. Weak lensing number counts covariance matrix for a ACDM 
model. The off-diagonal terms are small but not negligible only at low S/N 
ratios. The bottom panel shows the comparison between the diagonal com- 
ponents of the full covariance matrix derived numerically (shown in the top 
panel) and what would be expected for ideal number counts with Poisson 
statistic and statistically independent S/N bins. For high S/N the covariance 
is diagonal and Poissonian by construction. 



low S/N ratios up to S/N ~ 3, where we expect to have a small 
but not negligible covariance contribution since different S/N bins 
can be correlated. To account for the non-diagonal covariance el- 
ements, we realized a set of 1000 numerical realizations of weak 
lensing data, modeled with a Gaussian random field whose power 
spectrum, P(£), is discussed in Section l4!2l As expected the num- 
ber counts derived from the simulation agrees with the analytical 
predictions except for very low values of the S/N rati o (S/N < 0.2) 
where the analytic approximation fails as discussed bv lMaturi et al.l 
1(20100. In any case these very low S/N are completely negligible 
for any cosmological analysis. 

The total covariance matrix of weak lensing number counts is 
shown in Figure |4] As expected, it has a strong diagonal compo- 
nent and its off-diagonal terms are small but non negligible only at 
low S/N ratios. By treating the two components, LSS -1- noise and 
galaxy clusters as independent, we do not account for their mutual 
influence on the final number counts, but this effect is expected to 
be small because the signature of each cluster affects only few pix- 
els. In any case, this approach is conservative with respect to the 
final cosmological constraints because we neglect the impact of the 
detections with high S/N ratios, carriers of the non-Gaussianity sig- 
nal, on those with lower S/N ratios. 



6.2 Cosmological constraints on /nl 

In this Section we assume the expected performances of the weak 
lensing survey Euclid l iLaureiis..200 9). i.e. a field of view of 20, 000 
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Figure 5. Cosmological 1 - cr confidence levels as obtained for the local 
shape non-Gaussian model using weak lensing number counts, by assum- 
ing both Poissonian statistical independence (Poissonian), and by using the 
numerical simulations including the full covariance matrix (numerical). 



Figure 6. Confidence contours for the local shape non-Gaussian model, 
with different lower S/N cut-off limits. This shows the information lost by 
cutting the data at different S/N ratios. 
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Table 1. Dependence of the I - cr joint bounds on the background galaxy 
number density fjg (left) and on the filter scale radius, 6s (right). The fiducial 
model with erg = 0.809 and /nl = is adopted. 



deg^, an average background galaxy number density of = 40 
arcmin"^ with an intrinsic ellipticity nns of cr^^ =0.3, and a redshif t 
distribution similar to the one discussed bv lBeniamin et alj j2007l) . 
The confidence levels are evaluated with the Fisher matrix analysis 
discussed in Section 16.11 In Figure |5] we compare the confidence 
regions for the local shape non-Gaussian model obtained first by 
assuming purely Poissonian independent statistic, i.e. a data covari- 
ance having a pure diagonal component with Poissonian amplitude, 
and then by using the full covariance matrix obtained from the nu- 
merical simulations (see Section [Oi l. As can be seen, the 1 - cr 
confidence ellipse is slightly tilted in the latter case with respect 
to the former, implying a somewhat different direction of degener- 
acy between the two parameters under consideration. Overall, how- 
ever, the constraints on the non-Gaussianity level /nl are basically 
unchanged, while those on the normalization of the matter power 
spectrum erg are only very slightly loosened by using the simpler 
Poisson model. As a consequence, it is safe to state that the final 
differences between the two approaches are negligible and that for 
a first fast analysis the simplified approach can be used. 

Given that, we used the first approach to estimate the Cramer- 
Rao bounds on erg and /nl at 1 - cr level for different survey 
depths, i.e. different galaxies number densities, and for different fil- 
ter scales as shown in Table 1 . In order to save computational time. 



the full data covariance has been used only for the final results re- 
ported below. The results are evaluated for a survey sky coverage 
of 20, 000 deg^ and scale with the square root of the survey area. 
The constraints on both /ml and erg improve with increasing num- 
ber density of background sources. However only a marginal ben- 
efit is obtained by increasing the number density over 40 arcmin"^. 
In other words, in terms of joint constraints for the level of non- 
Gaussianity and amplitude of the matter power spectrum, it does 
not pay off to increase the depth of the survey beyond what is al- 
ready planned for Euclid. Concerning the scale radius of the op- 
timal filter, by increasing it we get a significant tightening of the 
constraints on /ml, while those on erg are almost insensitive to it. 

Now that we have an overall picture of the dependence on the 
survey characteristics, we used the numerical simulations includ- 
ing the full covariance matrix to obtain our final results. In Figure|6] 
we show the 1 - cr confidence levels in the erg - /^l plane ob- 
tained for the non-Gaussian model with local shape when different 
minimum S/N ratios are considered. It is in fact common practice 
for actual applications to use only those detections with S/N larger 
then S/N,ni„ = 3 - 5 in the attempt to isolate galaxy clusters from 
all other signal components. Here we show how the cut-off level 
affects the final results. The related loss of information happens be- 
cause the presence of non linear structures affects the statistic of 
weak lensing maps also at lower signal to noise ratios where they 
cannot be seen as clear single detections but where their presence is 
still visible. Moreover, although this effect is negligible for the /nl 
estimates, chance projection of the LSS dominates at these small 
S/N values, and still contains cosmological information. Figure |6] 
shows that adopting a conservative choice for the minimum S/N 
value can results in a loss of constraining power of a factor ~ 3 
on /nl and a factor of ~ 4 on erg. Thus, particular care have to be 
used when defining the detection selection criteria since data can be 
fully exploited only if lower S/N levels, e.g. S/Nn,in ~ 2, are used. 
This would be possible only if a deep statistical understanding of 
data is achieved together with a detailed modeling. A tentative to 
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confidence levels. 



go in this direction is provided by t he analytic predictio n for weak 
lensing number counts proposed bv lMaturi et alj falOl) but which 
still needs to improve some of the adopted approximations. 

Finally we show in Figures|7]the constraints for 68.3%, 95.4% 



and 99.7% confidence levels with S/N„ 



for the four dif- 



ferent primordial bispectrum shapes discussed in Section [2] We 
adopted the optimal filter with scale of 0s = 2'. As one could 
naively expect, while the constraints on erg are rather insensitive 
to the adopted shape of the primordial bispectrum, the constraints 
on /nl vary widely with it. Particularly, while the 1 - cr error A/nl 
is at the level of a few tens for the local and equilateral shapes, it 
grows up to ~ 100 - 200 for the enfolded and orthogonal shapes. 
These constraints are not competitive with what is expected by fu- 
ture CMB and galaxy clustering probes, and are only comparable 
with the expected pe rformance of the weak lensing power spectrum 
jpedeli & Moscardin i 2010) analysis. The main reason for this re- 
sults is that in this work we did not use the full redshift information 
which is otherwise included in the other works and wh ich could be 
inclu ded by the use of lensing tomography. Recently, iFedeli et"ail 
( I2OIO1) have shown that using the correlation function of weak 



lensing-selected clusters it is possible to push the bounds on /nl 
below what would be achieved without redshift information. In ad- 
dition, it is likely that a combination of these diverse probes can 
reduce the error on /nl substantially. We plan to investigate this 
synergy in a future work. 

The constraints on the amplitude of the matter power spec- 
trum, A cTg ~ 10~^ are them selves quite competitive. For instance, 
recentlv lWang et all j2010l) forecasted an error on erg of ~ 7 x 10"^ 
by using the shape and position of the Ba ryon Acoust i c Osc illation 
measured by Euclid. Despite the fact that lWang et al.l ( l201(t ) let all 
the cosmological parameters free to vary, while we limit our anal- 
ysis to /nl and erg, these numbers argue for our constraints on the 
amplitude of the matter power spectrum being sound. 

The degeneracy between the level of primordial non- 
Gaussianity and the amplitude of the matter power spectrum is also 
easy to understand, since increases in both /nl and erg bring an 
increment in the large scale matter power spectrum and the occur- 
rence of massive dark matter halos. The only exception to this is 
given by the orthogonal model, because in this case a positive /nl 
implies a decrement in both the abundance of cosmic structures and 
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Figure 8. The 1 - cr confidence levels in the erg - /nl plane for all the four 
shapes of the primordial bispectrum considered in this work, as labeled 

the large scale bias of dark matter halos. For further clarity, in Fig- 
ure d] we show the 1 - cr confidence levels in the erg - /^l plane 
for the four primordial bispectrum shapes considered in this work 
together. This Figure allows a more direct comparison between the 
different models, and highlights the differences in constraining the 
level of primordial non-Gaussianity between them. 



7 SUMMARY AND CONCLUSIONS 

The abundance of S/N peaks in cosmic shear maps contains cos- 
mological information through the mass function of large dark mat- 
ter clumps and the projection of the large scale matter distribution. 
We considered the effect of deviations from primordial Gaussianity 
on this particular weak lensing statistic. Analytic modeling of the 
galaxy cluster mass function, the large scale halo bias, and the fully 
non-linear power spectrum of dark matter was implemented, where 
the effect of primordial non-Gaussianity could be straightforwardly 
introduced. We predicted analytically and through simple numeri- 
cal simulations the expected weak lensing peaks number counts and 
covariance for different weak lensing filters. We then performed a 
Fisher matrix analysis in order to forecast the constraints on the 
level of primordial non-Gaussianity /nl and the amplitude of the 
matter power spectrum erg that could be expected by counting the 
cosmic shear peaks for different survey configurations, i.e. depth 
and field of view, in particular for the half-sky maps expected by 
the proposed ESA space mission Euclid. Our principal results can 
be summarized as follows. 

• Primordial non-Gaussianity affects mainly the high S/N part 
of the shear peak number counts, where the signal is dominated by 
the occurrence of real cluster-sized dark matter halos. For the local 
primordial bispectrum shape and according to the adopted weak 
lensing filter the number counts can be modified by ~ 3 - 4% at 
S/N > 10 if I/nlI = 50, and by up to ~ 10% for |/nlI = 100. For 
the enfolded and orthogonal bispectrum shapes, the modification 
reaches at most ~ 1 - 2% for |/nl| = 100. 



• Generically an increment in /nl corresponds to an increment 
in the number counts of cosmic shear peaks, and the other way 
round, compatibly with the effect of primordial non-Gaussianity 
on the cluster mass function. The exception to this is given by the 
orthogonal model, for which it is known that the skewness and 
scale-dependent bias correction are negative for positive /ml, and 
vice-versa. 

• We show that considering only cosmic shear peaks with S/N 
larger than 3 - 5, as it is costume in actual applications, leads to a 
loss of cosmological information because the LSS and, more im- 
portantly, the low mass clusters contributions are ignored. In par- 
ticular the latter is important in constraining /nl even if these struc- 
tures cannot be detected individually. Thus, data can be better ex- 
ploited if lower S/N levels, e.g. S/Nmi„ < 2, are used. For this 
achievement, a deep statistical understanding and a detailed data 
modeling are necessary. 

• In terms of joint constraints for /nl and erg only a marginal 
benefit is obtained by increasing the background galaxy number 
density above 40 arcmin"^, while it is more convenient to go for 
wider fields of view. By increasing the scale radius of the optimal 
filter we obtained a significant improvement for the /nl constraints, 
while those on erg are almost unaffected. 

• Counting cosmic shear peaks in future wide field optical/near- 
infrared surveys on the model of Euclid can constrain /ml to the 
level of a few tens for the local and orthogonal shapes, and to the 
level of ~ 100 - 200 for the enfolded and orthogonal shapes. Con- 
straints on CTg are at the level of ~ 10"^ and are instead rather in- 
sensitive to the assumed shape of the primordial bispectrum. 

The results presented in this work show that the abundance 
of cosmic shear peaks can be a powerful mean for constrain- 
ing cosmology. Forecasted errors on the level of primordial non- 
Gaussianity for a Euclid-like survey are competitive with similar 
constraints given by weak lensing tomography and the correlation 
function of weak lensing selected galaxy clusters. At the same time, 
counting cosmic shear peaks irrespective of their nature is much 
easier and less subject to systematics compared with other probes 
where the nature of each lensing peak, resulting by LSS or galaxy 
clusters, has to be determined. Our findings outline once more the 
great advance in understanding the physics of the primordial Uni- 
verse that is expected thanks to future wide field imaging siu-veys. 
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